Cosmological solutions of massive gravity on de Sitter 
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In the framework of the recently proposed models of massive gravity, defined with respect to a de 
Sitter reference metric, we obtain new homogeneous and isotropic solutions for arbitrary cosmological 
matter and arbitrary spatial curvature. These solutions can be classified into three branches. In the 
first two, the massive gravity terms behave like a cosmological constant. In the third branch, the 
massive gravity effects can be described by a time evolving effective fluid with rather remarkable 
features, including the property to behave as a cosmological constant at late time. 
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Long after the first attempt by Pauli and Fierz to give a mass to the graviton [l|, it has been realized, decades ago, 
that finding a healthy nonlinear massive extension of general relativity represents a formidable challenge because it 
requires to get rid of the so-called Bouldware-Deser ghost Q . Very recently, de Rham, Gabadadze and Tolley (dRGT) 
succeeded in constructing a massive theory of gravity that satisfies this criterion 0] , as later confirmed in 0| . Beyond 
■ its obvious theoretical interest, this achievement has a special significance in a context where most of the matter 
content of the Universe remains unknown and alternative explanations for dark energy and/or dark matter could 
reveal appealing. This explains why this recent model has attracted a lot of attention, especially for its cosmological 
consequences. In this respect, a surprising discovery was that dRGT massive gravity does not allow for spatially flat 
JL . homogeneous and isotropic solutions However, open cosmological solutions were obtained, with two branches of 
solutions in which the massive graviton terms lead to an effective cosmological constant [f| (other solutions relevant 
^ c~| for cosmology can be found in e.g. 04LTI)- 

In the present work, we start from a slightly modified version of the original dRGT massive gravity in which the (a 
| priori arbitrary) reference geometry is chosen to be de Sitter instead of Minkowski. A similar setting was explored very 
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recently in [12j and [13| . The de Sitter geometry possesses as many symmetries as the flat geometry but introduces a 
mass scale H c as additional parameter. In this setup, we have been able to find new cosmological solutions with flat, 
open or closed spatial geometry, for arbitrary cosmological matter. Our solutions can be classified in three branches, 
two of which being analogous to the open solutions of Ira , while the last branch exhibits a new and rich phenomenology. 
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Q 1 II. HOMOGENEOUS AND ISOTROPIC SOLUTIONS OF MASSIVE GRAVITY 

(N' 

We first present the theory of massive gravity introduced in Q , which can be described in terms of the usual four- 
dimensional metric g^ v and of four scalar fields <fi a (a — 0, ... ,3), called the Stiickclbcrg fields. Gravity is governed 
by the action 



ii? + m 2 (c,2 + ct 3 £ 3 + (x&Ca 



(2.1) 



where the first term is the familiar Einstcin-Hilbcrt Lagrangian (we set M p i = 1 in the following) and the three 
additional terms are specific functions of the metric <? M „ and of the four scalar fields </>°, via the tensor 



£"* = - Vfabg^d^d^" , (2.2) 

where f a b is called the reference, or fiducial, metric (the square root must be understood in a matricial sense). The 
explicit expressions for these additional terms in the Lagrangian are 

£ 2 = i([/C] 2 -[/C 2 ]) (2.3) 
^3 = i([/C] 3 -3[/C][/C 2 ]+2[/C 3 ]) (2.4) 
d = t^([£] 4 - 6[/C] 2 [/C 2 ] + 3[/C 2 ] 2 + 8[/C][/C 3 ] - 6[/C 4 ]) (2.5) 
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where the standard matrix notation is used (i.e. (/C 2 ) M „ = /C (T /C '„) and the brackets represent a trace. 

We now restrict our discussion to a FLRW (Fricdmann-Lemaitre-Robertson- Walker) geometry, of arbitrary spatial 
curvature, described by the metric 

ds 2 = g^dx^dx" = -N 2 (t)dt 2 + a 2 (t) - /ij (x)dx i dx j , (2.6) 

where the spatial metric 7^, written for example in terms of spherical coordinates, reads 

dr 2 

7« (x)dx*dxJ = — + r 2 (d9 2 + sin 2 6 def 2 ) (2.7) 

1 — kr z 

with k = 0, —1 or 1 for, respectively, flat, open or closed cosmologies. 

In the present work, we take for the reference metric f a b the de Sitter metric. As we will see, and in contrast with 
the Minkowski case, one can easily construct flat, open and closed cosmologies by starting from the appropriate slicing 
of de Sitter. Let us thus write the de Sitter metric in the form 

f ab dxj> a d(j) b = -dT 2 + b 2 (T) 7 y(X) dX l dX\ (2.8) 

where the functions b k (T) are defined by 

b (T) = e H " T , b-i(T)=H- 1 sinh{H c T), bx(T) = H~ l cosh(H c T) . (2.9) 

In the limit H c — > 0, one recovers the Minkowski metric in the flat and open cases: bo(T) = 1 and 6_i(T) = T, the 
latter case corresponding to the Milne metric for the flat geometry. 

We must now specify the Stiickclbcrg fields so that the cosmological symmetries are satisfied. One sees immediately 
that the choice 

4> a = T = f(t), <P=X l = x % (2.10) 
leads to a homogeneous and isotropic tensor, 

U» = fab d^ a d v 4> b = Diag {-/ 2 , b 2 (f(t)) 7ij } . (2.11) 
Denoting Sf the sign of /, the corresponding matrix /C, defined in (|2.2j) . is simply given by 1 

K,° = l-e f ^, K* s = (l - ^^Vi , /C l o = 0, !C° t = 0. 
Substituting in the Lagrangian of massive gravity, one gets 



(2.12) 



^9 = \f-g(£2 + CH3C3 + a 4 £ 4 ) 

= (a - h(f)) {N [a 2 (4a 3 + a 4 + 6) - a(5a 3 + 2a 4 + 3)6 fc (/) + (a 3 + a 4 )^(/)] 

-s } f [(3 + 3a 3 + a 4 )a 2 - (3a 3 + 2a 4 )ab k (f) + a 4 b k (f) 2 ] } . (2.13) 

The equation of motion for f(t) is obtained by varying this Lagrangian with respect to /: 

[(3 + 3a 3 + a 4 )a 2 - 2(1 + 2a 3 + a 4 )ab k (f) + (a 3 + a 4 )b 2 (f)] (± - e f b' k (f)j =0. (2.14) 
In general, there are several solutions for /. The first two solutions correspond to 



b k (f(t)) - X ± a®, X ± = 1 + 2Q3 + Q4±v ; i + Q3 + ^- a4 , (2.15) 

a 3 + a 4 

which exist only if the function b k is invertible. For a Minkowski reference metric f a b = Vab, one sees immediatly that 
there is no solution in the flat case since &o(£) = 1, whereas b-i(f) = f leads to two branches of solutions in the open 
case, in agreement with the conclusions of Q and [f|. 



1 We also assume / > in the case k = 



-1. 
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Let us now concentrate on the last branch defined by the condition 

e f b' k (f) = | • (2.16) 

It is non trivial only if b' k is an invertiblc function, which is not the case with a Minkowski reference metric, either in 
the flat or open cases. However, in our case, one can obtain an explicit solution for /(£) with the functions b k given 
in (|2.9[) . Before examining the flat case, let us stress that the solutions in this branch are necessarily accelerating as 
can be seen by taking the time derivative of (|2.16j) . which yields 

a = h'i{f) l/l > (N = l). (2.17) 

In the particular case k — 0, on which wc will focus in the following, one finds (assuming / > 0) 

where H denotes the usual Hubble parameter. 

III. FRIEDMANN EQUATIONS AND EFFECTIVE GRAVITATIONAL FLUID 

To obtain the Friedmann equations, one must add to C g the usual Einstcin-Hilbert term, which reads 

C E u = -^ + 3kNa, (3.1) 

as well as an arbitrary matter Lagrangian £ rn that describes ordinary cosmological matter. Variation of the total 
Lagrangian with respect to the lapse function N (which will be set to 1 in the following) then yields the first Friedmann 
equation 

3H 2 + 3^ = p m + p g , H=-, (3.2) 
ar a 

where p m denotes the ordinary matter energy density whilst p g corresponds to an effective energy density arising from 
the massive gravity action: 

o 
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p g = -§- (b k (f) - a) {(6 + ia 3 + a 4 ) a 2 - (3 + 5a 3 + 2a 4 )ab k {f) + (a 3 + a 4 )b 2 k (f)} . (3.3) 
The variation of the total action with respect to a(t) yields the second Friedmann equation in the form 



2H + 3H' + — = -P m - P g , (3.4) 

(7/ 



with the effective pressure 
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P g = ^ {(6 + 4a 3 + "4 - (3 + 3a 3 + a 4 )/) a 2 - 2 (3 + 3a 3 + a 4 - (1 + 2a 3 + a 4 )/) ab k (f) 

+ (l + 2a 3 + a 4 - (03 + a 4 )/) b 2 (f)} . (3.5) 
We now study the expressions of p g and P g for the three branches of solutions identified previously. 

1. First two branches 

Substituting the solution (|2.15|) . one finds that the massive gravity contribution behaves like a cosmological constant 
with 



1 + a 3 ± y/l + a 3 + d\ - a 4 ) ( 1 + a 2 - 2a 4 ± (1 + a 3 ) y/l + a 3 + a\ - a 4 

(a 3 + a 4 ) 2 



P g = -Pg = -TO 2 ^ : 7 \ — : 1 ■ (3.6) 



Note that the terms proportional to / in (|3.5p cancel because they are proportional to the combination that appears 
in the equation of motion for /. We recover exactly the result of [6|, even if the spatial curvature is no longer restricted 
to be negative. Remarkably, the result is independent of H c . 
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2. Third branch 



Let us now turn to the third branch where the effective gravitational fluid follows a much more sophisticated 
behaviour. For simplicity, we consider here only the flat case, but it is straightforward to extend the following analysis 
to the open and closed cases. Upon substituting the explicit solution (|2.18[) for / into (|3.3[) and (|3.5[) . one gets 



Pg 



H 



H H 

6 + 4a 3 + a 4 - (3 + 5a 3 + 2a 4 ) — + (a 3 + a 4 )-r^ 



(3.7) 



and 



P g = m 2 g 1 6 + 4a 3 + a 4 - (3 + 3a 3 + a 4 ) 
-(a 3 + a 4 )| 3 1 + ^) 




H 2 I H 
3 + 2- 



(3.8) 



It can be noticed that (|3.7[) coincides with the expression obtained by [13j in the special case of de Sitter cosmology, 
i.e. with a constant H , on a de Sitter reference metric, although the Hubble parameter is time-dependent in our 
case. One can check explicitly that the effective gravitational fluid, characterized by p g and P g , satisfies the usual 
conservation equation 



3H( Pg + Pg) =0. 



(3.9) 



The behaviour of the effective fluid described by the above energy density and pressure is quite peculiar. The energy 
density p g can be positive or negative and its sign can change during time evolution, when the ratio H J H c crosses 
specific values, which depend on the parameters a 3 and a 4 . For example, in the minimal model with a 3 = a 4 = 0, 
p g changes sign when H = H c or H = 2H C . 

Combining the second Friedmann equation with the first one yields the relation 



H 

3 + 3a 3 + a 4 - 2(1 + 2a 3 + a 4 )— + (a 3 + a 4 ) 

tin. 



H 2 



2R< 



H 

IT 2 



- (1 + w m )p m , 



(3.10) 



where w m 
rh 2 - 2H 2 



e P m /p m . One can thus identify a critical value for H when the term between brackets, of the form 
m 2 x(H/H c ) — 2H 2 , vanishes. Remarkably, this critical value coincides with the Higuchi bound derived 
in [13[ for de Sitter cosmology. Moreover, one can see that for the two other branches discussed earlier, the function 
x(H/H c ) vanishes, which seems to be related to the fact that the kinetic energy of the scalar mode around these 
solutions vanishes, as found in [l4j], since the kinetic term is proportional to m 2 (m 2 — 2H 2 ) according to [l3j]. If the 
Higuchi condition m 2 > 2iJ 2 is satisfied, the above relation implies that the Hubble parameter increases for matter 
satisfying the weak energy condition (i.e. w m > — 1 and p m > 0). 

Let us try to analyze the combined evolution of the effective gravitational fluid with ordinary matter, in the minimal 
model where a 3 = a 4 = for simplicity. It is convenient to introduce the dimcnsionlcss quantities 



Y 



H 



Pm 



A 



(3.11) 



Using the first Friedmann equation to express p m as a function of H, one finds that the above equation can be 
rewritten as a differential equation for Y and the full system can be written in the form 



Y' = 



3(1 + w m ) [(1 + A)F 2 - 3\Y + 2A] 
3A-2(1 + A)y 



= 3 [(1 + A)y 2 - 3\Y + 2A] 



(3.12) 



where a prime denotes a derivative with respect to the number of e-folds, i.e. Y = HY'. The second relation is simply 
a constraint between the value of the matter energy density and the Hubble parameter. In the following, we will only 
assume that the cosmological matter is characterized by r > and w m > — 1. It is then worth noting that the Higuchi 
condition m 2 — 2H 2 > corresponds to 



U = 3A - 2(1 + X)Y = -2(1 + A)(Y - Y H ) > 0, 



3A 



2(1 + A) 



(3.13) 
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In order to satisfy the Higuchi bound, one must therefore have Y < Yh if A > — 1, or Y > Yh if A < 
It is also useful to introduce the two roots of the numerator of the equation for Y, 



which are defined if A > 8 or A < 0. Rewriting the dynamical system (|3.12p in the form 

r =-|(l + ^ ^'y^-Yh^ ' r = '^ l + A )( y - Y +)( Y -Y-)>^ (3-15) 
it is easy to study its evolution, depending on the value of A: 

• A > 8 (which implies < F_ < Yh < Y + ): if H > 0, then Y < Y- and Y tends towards F_ asymptotically. By 
contrast, if the Higuchi bound is not satisfied, i.e. H < 0, one must have Y > Y + and Y decreases, converging 
asymptotically towards Y + . 

• < A < 8 (Y + and F_ are not defined) : the Higuchi bound is satisfied if Y < Yh initially, and Y increases to 
reach Yh in a finite time. By contrast, if Y > Yh initially, the Higuchi bound is not satisfied and Y decreases 
to reach Yh in a finite time. 

• — 1 < A < (Yh < 0): the Higuchi condition is never satisfied. The condition r > imposes Y > Y + and Y 
decreases toward Y + asymptotically. 

• A < — 1 (which implies < Y + < Yh < YL_): H > imposes Yh < Y < Y- initially and Y tends asymptotically 
towards Y~. If T~L < 0, one must have Y + < Y < Yh and Y decreases towards Y + . 

We thus find that in most cases (A < or A > 8), the effective gravitational energy density tends to a constant 
asymptotically, while the cosmological evolution approaches de Sitter, with a Hubble parameter that depends on A 
and is proportional to H c . When < A < 8, the system evolves towards a singularity at finite time. One can proceed 
similarly for general values 0:3 and 0:4 but the analysis is more complicated because the numerator and denominator 
of the equation for Y become respectively third-order and second-order polynomials in Y. 



IV. CONCLUSION 



In the present work, we have obtained spatially fiat (as well as open or closed) FLRW solutions with arbitrary 
cosmological matter in the context of ghost-free models of massive gravity, evading the no-go theorem of @ by 
adopting a de Sitter reference metric instead of Minkowski. The constraint equation for the Stiickelberg fields leads to 
three branches. In two branches, one finds that the effective gravitational fluid behaves like a cosmological constant, 
whose value, remarkably, is independent of H c and coincides exactly with the value obtained in Q for the specific case 
of open FLRW solutions with Minkowski as reference metric. By contrast, the third branch exhibits a much richer 
phenomenology, although expanding cosmological solutions are restricted to be accelerating. The massive gravity 
effects can be described by an effective fluid, which is in general time-dependent since its energy density depends on 
the physical Hubble parameter H (and its pressure on H as well). In the simplest case where 03 = 04 = we have 
investigated the cosmological evolution and found that the outcome is either a singularity at finite time or a de Sitter 
evolution, depending on the value of the ratio A = m 2 g / 'H%- 

To conclude, massive gravity on de Sitter leads to new solutions with surprising features. It would be worth exploring 
further these solutions, in particular by investigating more systematically the parameter space for the coefficients 013 
and CK4. It would also be interesting to study perturbations around these new solutions, by extending previous works 
on this topic (see e.g. fl3l. [L"H ) . 
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